Highly nonclassical photon statistics in parametric down conversion 
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We use photon counters to obtain the joint photon counting statistics from twin-beam non- 
degenerate parametric down conversion, and we demonstrate directly, and with no auxiliary as- 
sumptions, that these twin beams are nonclassical. 
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I. INTRODUCTION 

The quantum nature of light is exemplified by the 
fact that it can exhibit both undular (wavelike) and cor- 
puscular (particlelike) properties under different circum- 
stances. Using an atomic cascade to produce a pair 
of nearly simultaneous photons Grangier et al Q 
demonstrated that a heralded photon can either inter- 
fere with itself in a Mach-Zehnder interferometer or be 
shown to either wholly transmit or reflect by a beam 
splitter, which demonstrates complementarity for light. 
More recently the continuous trade-off of undular and 
corpuscular features in a quantum non-demolition mea- 
surement was demonstrated using a non-deterministic 
linear optical gate p|. Technically, nonclassical light 
corresponds to field states whose Glauber-Sudarshan P- 
representation is singular 0, [|| ; practically nonclassical- 
ity can be established by demonstrating photon anti- 
bunching or using optical homodyne tomography 
to demonstrate that the Wigner function (which is a 
Gaussian convolution of the P-representation) exhibits 
squeezing |lOj| or is negative over some region 

Recently Waks, Diamanti, Sanders, Bartlett, and Ya- 
mamoto (WDSBY) ^| exploited cutting edge photon 
counting technology (using the Visible Light Photon 
Counter [l^) to demonstrate that a single beam emit- 
ted by parametric down conversion (PDC) is nonclassi- 
cal by measuring its photon statistics and avoiding aux- 
iliary assumptions such as detector efficiency. Here we 
reinterpret the WDSBY data using Klyshko's powerful 
criterion for direct detection of nonclassicality via pho- 
ton statistics 14], and we demonstrate nonclassicality of 
twin beam PDC by using two photon counters and Lee's 
extension of Klyshko's criterion to two beams [l5| . 

Although the simultaneity of photon pair produc- 
tion by twin-beam, or nondegenerate, PDC is well 
known [l6| , leading to its applications as a heralded pho- 
ton source an d f° r tests of Bell's inequalities |18|. 
the joint photon statistics of the twin beams have not 
previously been directly measured. These statistics are 
necessary for a direct confirmation of the nonclassicality 
of twin-beam PDC, which is one of the most important 
tools in quantum optics as a source of correlated or en- 
tangled pairs of photons |l9l | and, more recently, as a 



source of entangled four-tuples of photons 20] . We have 
performed this direct photon counting experiment and 
obtained joint photon counting data for the twin beams 
and thereby showed that these twin beams are truly non- 
classical. 



II. SINGLE BEAM NON-CLASSICAL 
STATISTICS 

Theoretical predictions suggest that the beams should 
be nonclassical, but theories generally treat the twin- 
beam PDC output as a two-mode squeezed vacuum 
state 
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with rj G C and the argument of rj dependent on the phase 
of the pump laser field for the PDC. For degenerate PDC, 
the state is replaced by 
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and the photon number distribution is 

^ = (1-|X| 2 )W 2 ". (3) 

This state corresponds to a joint photon number distri- 
bution for the two beams, which is given by 

Pn u n 2 =(l-\rf)\ri\ 2ni Sn u n 2 . (4) 

Klyshko considered a state of the type emitted by an ideal 
degenerate PDC, and Lee's analysis of photon counting 
correlations explicitly assumed two-mode states and in- 
cluded a treatment of the two-mode squeezed vacuum (Q . 

The reality of PDC and photon counting is somewhat 
more complicated. The twin beams (known as signal and 
idler fields) are each multimode so the theoretical treat- 
ments of these outputs are not immediately applicable. 
In fact there are so many modes for each of the signal and 
idler beams that we observe a photon distribution over 
these modes that is Poissonian. The Poisson distribution 
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originates from the fact that there so many signal and 
idler modes that the average number of pairs per mode 
is much less than one. This means that each mode is 
predominantly a vacuum state with a small one-photon 
contribution and a negligible multi-photon contribution. 
Thus, the generated signal is created by a sum of inde- 
pendent spontaneous emitters (one for each mode), and, 
as each pair is created independently, the pair creation 
statistics must be a Poisson distribution. 

This discrepancy in photon statistics between the usual 
theoretical treatment of PDC output [2l| and experimen- 
tal reality makes the direct measurement of joint photon 
counting statistics and verification of nonclassicality even 
more interesting and important. Properties of twin-beam 
PDC have been explored in an elegant experiment by 
Haderka et al (2^ , but their approach is different in that 
they do not use two distinct photon counters on each 
beam. Thus, inferences of nonclassicality rely on auxil- 
iary assumptions that are not required by our approach 
such as taking into account the losses during transmis- 
sion, quantum efficiency and internal noise of the camera 
and noise due to other light sources to obtain the joint 
signal- idler photon number distribution p2^|. 

To begin, we re-analyze the WDSBY data for single 
beam PDC. WDSBY [lj] conclusively demonstrated non- 
classicality using their data for one-, two-, and three- 
photon measured counts {p„;n = 1,2,3}, where we use 
the symbol p for the ideal probability and p for the mea- 
sured probability. We will now compare their results 
to Klyshko's criterion We introduce Klyshko's cri- 
terion for two reasons: (i) Klyshko's criterion demon- 
strates that the WDSBY data violates classical bounds 
beyond the three-photon case studied by WDSBY, and 
(ii) Klyshko's criterion is a basis for Lee's criterion 0] 
which we use for the analysis of the twin beams. 

WDSBY used the fact that nonclassicality of light, 
defined by the singularity of the Glauber-Sudarshan P- 
representation holds if and only if the photocount dis- 
tribution for the beam cannot be expressed as a sum 
or integral of Poisson distributions (which corresponds 
to classical count distributions). Specifically the photon 
count statistics {p n ; n = 1, 2, 3} were irreconcilable with 
classical photon statistics; thus photon counting provided 
a direct means for establishing that a single-beam light 
source is nonclassical. 

To demonstrate nonclassical photon statistics, 
WDSBY pumped a Type-I phase-matched BBO crystal 
set up for collinear degenerate amplification with 20ns 
pulses of the fourth harmonic (266nm) of a Q-switched 
Nd:YAG laser. In this configuration, the down-converted 
photons had half the energy of the pump (532nm) and 
travelled in the same direction. The pump was removed 
by a prism, while the down-conversion was focused by 
a 250mm lens onto the VLPC detector. The detector 
output was amplified and then sent to a gated boxcar 
integrator, which was triggered by the laser. This 
configuration was used to measure the photon number 
distribution of the detected field. WDSBY achieved 



a forty standard deviation violation of classicality, 
according to the criterion 
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0.379 (5) 



without any auxiliary assumptions about the detector or 
field. 

In fact Klyshko's criterion preceded the WSBY crite- 
rion and was the first proposal for a direct test of non- 
classicality according to 'local properties' of the photon 
number distribution (PND) 



al pr 



K„ 



(n + l) P "- lP " +1 <1, n= 1,2,3, .... (6) 



By the replacement p t— > p, the criterion applies to 
measured photon statistics rather than the ideal photon 
counting distribution. Lee refers to such tests of non- 
classicality according to local properties of the PND as 
Type II, as opposed to the traditional version (which Lee 
calls Type I) that employs inequalities on moments (such 
as the Mandel Q parameter |23j). 

If criterion © is satisfied for any n, then the field is 
necessarily nonclassical. To compare with T, consider 
n = 2 for Klyshko's criterion applied to measured data: 
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Rearranging the terms in Eq. Q yields the nonclassical- 
ity criterion 
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An alternative criterion emerges from the nonclassicality 
criterion JSJ), which yields 
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Thus, combining Klyshko's result for K2 and the 
WDSBY result for T, we obtain the general condition, 
based on measured results {pi, p2, ps} for nonclassical- 
ity being 



p 2 > min-j \/^Pip3,^y^(pi - v'::> 



(10) 



This condition is based on the best case by combining 
the WDSBY T criterion (JSJ with Klyshko's criterion for 
n = 2. 

Eq. I|10|) can be further simplified: the two quantities 
in this equation coincide in the symmetric case that pi = 
P3: then p2 > \Jij1p\. Now suppose that p3 = Kp\. 
We can thus rewrite condition (jl0|) as 
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FIG. 1: Experimental results for Klyshko's figure of merit K 
vs photon number n for a single-mode field. 



Except for k = 0, 1 the first term on the left-hand side of 
Eq. Ijlll) is always smaller than the second term, and the 
two terms are equal when k = 0, 1. Thus Klyshko's condi- 
tion is stronger than the WDSBY criterion, and another 
advantage of Klyshko's general criterion over WDSBY's 
criterion is that the latter only applies in the measured 
photon number distribution localized around n = 2 
whereas Klyshko's criterion explores nonclassicality for 
local regions of the photon number distribution around 
any n > 1. 

In Fig.^ we present a plot of K n vs n for the WDSBY 
data. Klyshko's criterion is violated for even numbers of 
photons between n = 2 and n = 8. The limit of our mea- 
surement capability is nine photons, and all odd-number 
detections fail to violate Klyshko's criterion as expected 
from joint distribution fllf. The reason why even- number 
counts violate the classical limit and odd-number counts 
do not is that the ideal photon number distribution is 
non-zero only for even- number counts. Experimentally 
we detect odd-number counts due to finite detector quan- 
tum efficiency, but this does not lead to violation of the 
classical limit. Fig. 1 is especially important in demon- 
strating that nonclassicality of the beam is evident for 
all even photon numbers up to the limits imposed by our 
photon counter capabilities. 



III. TWIN BEAM NON-CLASSICAL 
STATISTICS 

Now let us consider the twin-beam PDC and its joint 
photon counting statistics. We will use Lee's generaliza- 
tion of Klyshko's Type II criterion to establish nonclas- 
sicality of the twin-beam PDC output directly from the 
joint photon counting data with no adjustments made for 
photodetection efficiency. We will need to adapt Lee's 
criterion to accommodate the multi-mode aspect of the 
PDC output. 

Lee introduces the non-classical criterion for two 



modes of a light field given by 
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with 7ii,ri2 € {1,2,3,...}. The field is nonclassical if 
Rn!n 2 satisfies inequality (|12fl for any n\ and n 2 . Al- 
though Lee's analysis supposes a two-mode field, this re- 
sult is equally valid for a multimode field and the use of 
two photon counters, with some modes directed to one 
detector and other modes directed to the other detector. 

Mode mismatches and lost modes are partially respon- 
sible for photon counter inefficiencies. As we do not em- 
ploy auxiliary assumptions, the onus is on us to enhance 
efficiency to ensure that inequality (|12f> is satisfied, and 
not to resort to modeling these losses and adjusting the 
photon statistics according to these assumptions. 

We define n\ as the number of signal photons and n 2 
as the number of idler photons. It is convenient to deter- 
mine the conditional probabilities for n 2 photons given 
a count of n\ at the other detector, denoted by p n2 \m- 
The conditional and joint number counting probabilities 
are related by the formula 



Pn 1 ,n 2 — Pn 2 \n 1 Pn 1 ' 



(13) 



As explained in the previous section, the counting statis- 
tics for each detector is Poissonian distribution with 



p n = exp(-n)n n /nl 



(14) 



where n is the mean photon number summed over all the 
signal modes. 

The choice of the poisson distribution is justified by 
the fact that there are many more signal and idler modes 
than photons distributed amongst these modes [22| . 
Thus, each mode is predominantly a vacuum state with 
a small one-photon contribution and a negligible multi- 
photon contribution. The photon distribution across all 
modes of each beam is therefore binomial, which, in the 
limit of a large number of modes and fixed number of 
photons, is Poissonian. This Poissonian distribution con- 
trasts with the 'thermal' photon statistics for an out- 
put mode of ideal nondegenerate PDC, which follows 
from JTJ. 

Lee's expression (|12f) can be revised in terms of 
conditional probabilities. Using the simple relation 
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/ Pn-T° n = n/n, we rewrite the criterion as 



n 2 p„ 2 _i|„ 1+ i + n\{ni + l)p 

n 2 + l\ni- 

2nn 2 p n2 , ni 



< 1. 



( 15 ) 

In our experiments, the mean photon number collected 
over all the modes is n = 1. 

The experimental setup for generating twin photon 
beams for measuring the Lee non-classical criterion is 
shown in Fig. El A 266nm pump source is generated from 
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FIG. 2: Experimental setup for measuring violation of Lee's 
inequality. 

the fourth harmonic of a Q-switched Nd:YAG laser. The 
pump pulses have a duration of 20ns, and a repetition 
rate of 45kHz. A dispersion prism separates the fourth 
harmonic from the residual second harmonic, which is 
used to illuminate a high speed photodiode to generate 
a triggering signal. The fourth harmonic pumps a BBO 
crystal, which is set for non-collinear degenerate phase- 
matching. 

In this condition, the signal and idler waves are both 
532nm in wavelength and have a divergence angle of 1 
degree from the pump. The pump is loosely focused be- 
fore the BBO crystal to achieve a minimum waist at the 
collection lens. This results in a sharper two-photon im- 
age which enhances the collection efficiency. The pump 
power is set to 20/^TV. Using the count rate on the detec- 
tors and known values of the detector quantum efficiency, 
the average pair creation rate at this pump power is mea- 
sured to be one pair per pump laser pulse. 

Two VLPC detectors are used in this experiment. 
Each detector is held in a separate helium bath cryostat 
and cooled down to 6-7K, which is the optimum operat- 
ing temperature. The VLPC is sensitive to photons with 
wavelengths of up to 30/mi, so it must be shielded from 
room temperature thermal radiation. This is achieved by 
encasing the detector in a copper shield, which is cooled 
down to 6K. Acrylic windows at the front of the copper 
shield are used as infra red filters. These windows are 
highly transparent at visible wavelengths and simultane- 
ously nearly opaque at 2-30/mi wavelengths. 

VLPC 1 is used as the triggering detector which de- 
tects the number of photons generated in the signal arm 
on a given laser pulse. The output of VLPC 1 is amplified 
by an integrating amplifier, which generates an electrical 
pulse whose height is proportional to the number of emit- 
ted electrons. The height of the pulse is discriminated by 
a single channel analyzer (SCA). A logical AND is per- 
formed between the output of the SCA and the output of 
the photodiode to reject all detection events which occur 
outside of the pulse duration. 

The output of VLPC 2, which measures photons in 
the idler arm, is amplified and connected to the signal 
input of a boxcar integrator. On each pulse from the 
SCA, the output of VLPC 2 is integrated over a 2/is win- 
dow, which is sufficiently large to encompass the entire 



electrical pulse (determined by the bandwidth of subse- 



TABLE I: Experimental results for Lee's figure of merit R vs 
photon number n for twin beams 
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0.69 ± 0.023 


0.27 ± 0.007 


0.47 ± 0.012 


1.37 ± 0.04 
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2.23 ± 0.08 


0.70 ± 0.022 


0.33 ± 0.008 


0.47 ± 0.012 



quent amplifiers). The pulse area is proportional to the 
number of detected photons. By measuring the pulse 
area histogram of the VLPC, we can therefore measure 
Pn 2 \n 1 1 where n\ and ni are the number of photons in the 
signal and idler arm respectively. 

The experimental results are given in TableQ] Since we 
can only measure ri\ from one to four photons, we only 
obtain expressions for the Lee bound for n\ — 2, 3. In 
the ideal case where wc have perfect detection efficiency 
and no dark counts, we would have R ni n 2 = when 
ni = 7i2- When rix ^ n-i the Lee bound would not be 
well defined because p ni ,n 2 = ; causing a divergence. 
In the presence of detection losses, we would still expect 
the best violations when m = Ti2, with worse results in 
the off-diagonal term. As can be seen in Tabled this 
is indeed the case. When m — ri2 G {2,3}, we obtain 
extremely good violations of the Lee criterion, whereas 
off-diagonal terms yield worse violations or no violations 
at all. 



IV. CONCLUSIONS 

In summary we have demonstrated directly that twin- 
beam PDC produces non-classical light. As PDC is one 
of the most important tools for quantum optics, it seems 
surprising that the nonclassical nature of PDC output 
has not been directly verified before. One reason for 
not having previously establishing nonclassicality of twin- 
beam PDC is the requirement for sophisticated, modern 
photon counters and correlations of their data. Another 
reason for the novelty of our results is that only recently 
has direct testing of nonclassicality via the local proper- 
ties of the measured photon count distribution been un- 
derstood 

111 13 113 Our method of directly measuring 
photon counts for twin beams applies to multiple-beam 
fields and testing the local properties of the measured 
photon count distribution provides a valuable, practical 
means for establishing nonclassicality of light, especially 
in cases where the photons are not anti-bunched. 
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